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ABSTRACT 

We use a combination of high-resolution gasdynamics simulations of high-redshift dwarf galaxies and dissi- 
pationless simulations of a Milky Way sized halo to estimate the expected abundance and spatial distribution of 
the dwarf satellite galaxies that formed most of their stars around z ~ 8 and evolved only little since then. Such 
galaxies can be considered as fossils of the reionization era, and studying their properties could provide a direct 
window into the early, pre-reionization stages of galaxy formation. We show that ~ 5- 15% of the objects 
existing at z ~ 8 do indeed survive until the present in the Milky Way like environment without significant evo- 
lution. This implies that it is plausible that the fossil dwarf galaxies do exist in the Local Group. Because such 
galaxies form their stellar systems early during the period of active merging and accretion, they should have 
spheroidal morphology regardless of their current distance from the host galaxy. Their observed counterparts 
should therefore be identified among the dwarf spheroidal galaxies. We show that both the expected luminosity 
function and spatial distribution of dark matter halos which are likely to host fossil galaxies agree reasonably 
well with the observed distributions of the luminous (Ly > 10 6 L©) Local Group fossil candidates near the host 
galaxy (d < 200 kpc). However, the predicted abundance is substantially larger (by a factor of 2-3) for fainter 
galaxies (Ly < 10 6 L Q ) at larger distances (d > 300 kpc). We discuss several possible explanations for this 
discrepancy. 

Subject headings: cosmology: theory - galaxies: dwarf - galaxies: evolution - galaxies: formation - stars: 
formation - galaxies: halos - methods: numerical 



1. INTRODUCTION 

Abundance and spatial distribution of the faintest dwarf 
galaxies provide information and constraints on the galaxy 
formation processes and may give important cl ues to the na- 
ture o f dark matter. Cosmological simulations ( Rlvp inet alJ 
19991 iMoore et alJ 11999ft and semi-analytic mod els of 
galaxy formation t Kauffmann et aH [1993: Bullock et al. 
l2000nSomervillel2002UBenson et alJl20 02) predict the num- 
ber of gravitationally bound, dark matter clumps (or subhalos) 
around the Milky Way type galaxies larger, by up to an order 
of magnitude, than the observed number of dwarf galaxies in 
the Local Group. In addition, the predicted spatial distribu- 
tion of the dark matter subhalos around their hosts is more 
extended than the observed distribution of the Local Group 
satellites (Tavloret al. 2004; Rravtsov et al. 2004), indicating 
that strong spatial bias in galaxy formation should have been 
present if the CDM scenario is correct. 

Formation of galaxies in small-mass, dwarf-sized halos 
should be affected by heating from the ultra-violet (UV) radi- 
ation and internal feedback from stars. UV heating, in particu- 
lar, is expected to be highly important after the epoch of reion- 
ization when the entire volume of the Universe is affected by 
the cosmic UV background. The UV heating leads to dra- 
matic increase of the characteristic Jeans or filtering mass, 
M F , of intergalactic gas JGnedin & Huilll998tlGnedinll2000h . 
The dark matter halos with masses M < Mp (or circular ve- 
locities V c < Vf) cannot accrete new gas and therefore do not 
have supply of fresh gas to fuel the continuing star formation 
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(e.g., Efstathiou |[T99l iThoul & Weinberdfi996t lOuinn et al] 
| 1996llNavarro & Steinmetzll997llKitavama & Ikeuchl200Ct 
ICmedrnl l2000t iDiikstra et al] 12004^ In addition, the UV 
backg round can evapo r ate the existing gas in small halos 
( Bark ana & Loer3fi999T: IShaviv & Dekelll2003t IShapiro et alJ 
120041) . which would also reduce the star formation rate. 

In the presence of these galaxy formation suppression pro- 
cesses, the smallest dwarf galaxies can form: (1) if their 
host halos have formed sufficiently early, befor e the reioniza- 
tion e poch and UV heating became significant ( B ullock et al] 
2000) and/or (2) if the host halos become sufficiently massive 
after reionizati on such that UV heating effects are not signif- 
icant ( Rravtsov et al.l l2004l) . It is of course possible that a 
halo forms early and forms stars before reionization and then 
stays above filtering mass for an extended period of time after 
reionization due to continuing merging and accretion. Such 
gala xies would exhibit c ontinuous, albeit perhaps episodic 
(see [Rravtsov et al. 2004), star formation over an extended 
period of time. 

It is also possible that some galaxies form stars before 
reionization and then not evolve much in terms of their to- 
tal and stellar mass, with their stellar populations evolv- 
ing mostly passively. Most of such galaxies do not sur- 
vive until z = but are disrupted via merging and tidal 
heating, contributing to the formation of the bulge, stel- 
lar halo, and halo populatio n of globular clusters of the 
host galaxies (e.g.. Iwhite & Springe! I2000t iBullock et ap 
1 200 11 iRravtsov & Gnedinl l2005t iBuUock & Johnstonl l2005t 
Robertson et al. 2005; Moore et al. 2005^. Howe ver, a frac- 
tion of them can survive. Ricotti & Gnedin (2005) found that 
many of the dwarf galaxies identified in the high-resolution 
gas dynamic simulation at z ~ 8 do indeed resemble a sub- 
set of the observed old dwarf spheroidal galaxies in the Local 
Group in their kinematic, structural, and chemical properties. 

Such galaxies would then represent fossils of the reioniza- 
tion era, and studying their properties could provide a direct 
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window into the early, pre-reionization stages of galaxy for- 
mation. As such, it would be extremely interesting to identify 
and study this population in the Local Group. In this paper, 
we present theoretical estimates for the abundance and spa- 
tial distribution of the fossil galaxies using a combination of 
sophisticated galaxy formation simulations at high redshifts, 
which include non-equilibrium chemistry and self-consistent, 
three-dimensional radiative transfer and high-resolution dis- 
sipationless simulations of a Milky Way-sized halo, which 
allows us to track small-mass dwarf dark matter halos from 
the pre-reionization epoch (z ~ 10) to their present-day Milky 
Way-like environment. 

2. METHOD 

2.1. Simulations 

In this paper we combine results from several simulations 
to connect models of high-redshift dwarf galaxies with the 
present-day Local Group. The simulation we use for mod- 
eling high-redshift dw arf galaxies is described in detail in 
Ricottiet alJ J200 2albl) a s run "256Llp3" and was used in 
Ricotti & Gnedin (2005) as the main basis for their conclu- 
sions. 

The simulation includes dark matter, gas, stars, the feed- 
back of star formation on the gas chemistry and cooling, 
and the spatially-inhomogeneous and time dependent radia- 
tive transfer of ionizing and Lyman- Werner band photons - the 
latter being crucial for realistic modeling early star formation. 
The run followed 256 3 dark matter particles, an equal num- 
ber of baryonic cells, and more than 700,000 stellar particles. 
The mass resolution of the simulation is 900 M Q in baryons, 
and real comoving spatial resolution (twice the Plummer soft- 
ening length) is 150/i _1 pc (which corresponds to a physical 
scale of 24 pc at z = 8.3) in a computational box with of 
1/T 1 Mpc on a side. Hereafter, we refer to this simulation 
as LI. 

Clearly, the \h Mpc box is too small to be a representative 
volume of the universe at any moment in time. The box is ap- 
propriate for the purposes of this study because properties of 
the small-mass dwarf galaxies formed in this simulation are 
determined primarily by the physics of gravitational collapse, 
gas cooling, star formation and feedback and should not be af- 
fected by the box size. However, these galaxies would not be 
able to reionize even such a small volume of space. Therefore, 
iRicotti & Gnedinl (|2p05 ) introduced an additional source of 
ionizing radiation within the computational box, correspond- 
ing to a star-forming galaxy with constant star formation rate 
of lM©/yr. The source was switched on at z = 9.0, and by 
z = 8.3 the entire simulation box was completely ionized. At 
that time, the computational volume contains about two dozen 
galaxies with luminosities between 10 5 and 1O 7 L and circu- 
lar velocities < 30km/s. These galaxies l ost their d iffuse IGM 
due to the UV feedback of reionizatio n JGned inl l2000l) . and 
stopped forming stars. Ricotti & Gnedin (2005) argue that 
these galaxies can be identified with the population of fossil 
galaxies in the Local Group dwarfs. 

The specific value for the redshift at which 1/T 1 Mpc box 
is reionized (z = 9) is, of course, somewhat arbitrary, but it is 
important to note that it should not be identified with the red- 
shift of reionization of the entire universe. It is likely that the 
region around the Local Group progenitor is reionized consid- 
erably earlier than the universe as a whole, because progen- 
itors of the Milky Way-sized objects can easily create their 
own H II regions of the order of 1/T 1 Mpc by z ~ 9 or even 



earlier . 

A 1/T 1 Mpc box size is also too small to be representative 
of the Local Group environment. In the concordance LCDM 
cosmology 5 adopted throughout this paper, this box contains 
only 1.2 x 10 u Mq, whereas the mass o f the Local Group i s 
estimated to be about 3 x 10 12 M Q (e.g.. lKlvpin et aU l2002). 
In order to make a connection between z = 8 universe and the 
Milky Way-like environment at z = 0, we use a dissipation- 
less (dark matter only) simulation with 20/T 1 Mpc box (here- 
after L20), which was run using the Adaptive Refinement Tree 
JV-body code (ART, iKravtsov et al Jl 1 9971 iKravtsovl 1 9 99 ) to 
follow the evolution of the Milky Way-sized halo with high 
resolution. A Lagrangian region corresponding to a sphere of 
5/? v i r of a Milky Way-sized halo at z = in a low-resolution 
simulation was resimulated with higher resolution using mul- 
tiple mass r esolution tec hnique (this simulation is Box20 run 
presented in lPrada et alJl200 5f). In the high-resolution region 
the mass of the dark matter particles is m p = 6. 1 x 10 5 /z _1 M , 
corresponding to effective 1024 3 particles in the box, at the 
initial redshift of the simulation (zi = 70). The high mass res- 
olution region was surrounded by layers of particles of in- 
creasing mass with a total of four particle species. The sim- 
ulation starts with a uniform 256 3 grid covering the entire 
computational box. Higher force resolution is achieved in 
the regions corresponding to collapsing structures by recur- 
sive refining of all such regions using an adaptive refinement 
algorithm. Only regions containing highest resolution parti- 
cles were adaptively refined. The maximum of nine refine- 
ment levels was reached in the simulation corresponding to 
the peak formal spatial resolution of 150/T 1 comoving par- 
sec. The Milky Way-sized host halo has the virial mass of 
1.4 x 10 12 /! _1 Mq (or 2.3 million particles within the virial ra- 
dius) and virial radius of 230/T 1 kpc at z = 0. 

2.2. Simulation and Data Samples 

Before we compare results of the cosmological simulation 
with observational data, we consider the limits of its numeri- 
cal resolution. In this paper we choose to characterize dark 
matter halos by their maximum circular velocity V max , be- 
cause it is well-defined while definition of the halo mass is 
rather arbitrary. For the low-mass halos that may host dwarf 
galaxies, the maximum circular velocity function (the num- 
ber of halo s per unit V max ) is expected to be close to a power- 
law (e.g.. lKlvpin et aljl!99 9). Thus, strong deviation from a 
power-law at small halo masses indicates resolution limit of 
the simulation (see, e.g. JReed et alJ2005l) . 

Figure^iHustrates this point. It shows the cumulative max- 
imum velocity function (i.e., the number of dark matter halos 
with y max above a given value) at z = 8 and z = for the L20 
run. The strong deviation from the power-law at z = 8 can 
be seen at V max < 13km/s, indicating the resolution limit of 
the L20 simulation at this epoch. At z = 0, the issue is more 
complicated because resolution depends on the distance to the 
Milky Way-sized halos halo. For example, an isolated halo 
will exist even if its internal structure is poorly resolved, but 
it may be prematurely disrupted when it falls into a parent 
halo and experiences tidal heating. Three lines in the lower 
panel of Fig.^do indeed show that while the L20 simulation 
is sufficient to resolve dark matter halo down to about 8 km/s 

4 This indeed happens in a gasdynamics simulation of a Milky Way-sized 
progenitor with radiative transfer simulation analyzed by the authors (the sim- 
ulation is described briefly in Tasitsiomi 200a). 

5 (n m , U A , h, <x 8 ) = (0.3, 0.7, 0.7, 0.9). 
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FIG. 1. — Cumulative maximum circular velocity functions for the L20 
run at z = 8 (top panel) and z = (bottom panel). Velocity functions at z = 
are shown for three subsamples: halos within 1 Mpc of the parent halo (light 
gray line), 300 kpc (dark gray line), and 100 kpc (black line). Tilted thin 
dotted lines show power-law N(> V m ax) on V,^ ax , while the vertical dotted 
lines mark characteristic values of V raax corresponding to the resolution limit 
of the simulation, as discussed in the text. The velocity functions in both 
panels are normalized to unity at V max = 20 km/s. 



within 1 Mpc of the parent halo 6 , it is complete only down to 
17 km/s for halos within 100 kpc of the Milky Way-sized halo. 

We can simply correct for the incompleteness by a velocity 
dependent factor that is required to restore the power-law form 
of iV(> Vmax) function (in other words, the correction factor 
is the ratio of the power-law and the actual velocity function 
measured in simulations). To illustrate the effect of incom- 
pleteness we will use both the uncorrected and the corrected 
samples, which should bracket the correct result. As we show 
below, our results are not highly sensitive to this correction 
even for the d < 100 km/s sample, for which the correction is 
the largest. Note also that the correction only becomes large 
at y max 15 km/s, while none of the observed dwarf galax- 
ies in the Local Gro up are likely to have V max significantly 
below 15 km/s (e.g.. | Stoehr et aL]|2002{: |Zentner & Bullock 
120031 iKazantzidis et alJ2004tlMoore et alJ2 005 ) 

Figure |2] shows the sample of dark matter halos from the 
L20 run we use in this paper. Three shades of gray show 
three subsamples according to the distance from the parent 
halo at z = 0. As we mentioned above, we only select halos 
that have V max > 13 km/s at z = 8, so our sample is complete 
at high redshift. For these halos, V max at z = reaches down 
to Vmax ~ 8 km/s, so our sample of halos with d < 1 Mpc is 
effectively complete, while subsamples with smaller limits on 
d are incomplete to various degrees. 

Finally, the LI simulation allows us to relate the stellar lu- 
minosities of model dwarf galaxies to the maximum circular 
velocities of their dark matter halos at z = 8. The latter are 
measured with the same halo finder algorithm used to iden- 
tify halos in the L20 run. Figure [3] shows the distribution of 
these two properties for the model dwarf galaxies from the 
LI run. We marked the completeness limit of the L20 run at 

6 The improvement in resoltion limit, from 13 km/s at z, = 8 to 8 km/s at 
z = 0, is likely due to the fact that in a adaptive mesh refinement simulation 
we use, the resolution becomes somewhat better with time. 
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FIG. 2. — Distribution of maximum circular velocities for dark matter halos 
from run L20 at z = 8 versus z = 0. Halos at z = are distinguished by their 
distance to the parent halo: halos within 1 Mpc (large light gray squares), 
within 300 kpc (intermediate darker gray squares), and within 100 kpc (small 
black symbols). Only halos with V max (z = 8) > 13 km/s are shown. 



z = 8 (Vmax = 13 km/s) with the thin dotted line, and throughout 
this paper we only use dwarf galaxies above this limit. For- 
tunately, all model dwarfs with Ly > 10 5 L Q fall within this 
range. 

The objects shown in Figures |2 and [3J thus, constitute 
the sample of the simulated galaxies that we use in this pa- 
per. Our main task now is to compare these model galax- 
ies to the observed dwa rf galaxies in the Local Group that 
iRicotti & Gnedinl (|2005) identified as reionization fossils. 
The observed sample includes Draco, Phoenix, Sculptor, Sex- 
tans, Tucana, and Ursa Minor for the Milky Way subgroup, 
and And I-ELV.VIJX, Antlia, Cassiopea, Cetus, EGB0427, 
and SagDIG for the Andromeda subgroup. These galaxies 
are selected because they are dim (V-band luminosity below 
1O 7 L ) and because they formed at least 70% of their stars 
in one single burst at an early epoch (which is identified with 
the reionization epoch). These galaxies therefore represent 
reasonable candidates for the reionization fossils. This iden- 
tification is, of course, somewhat uncertain because the ab- 
solute ages of 10- 12Gyr old stars are not measured with 
sufficient accuracy. Nevertheless, we will assume in the re- 
mainder of this paper that the observational samples are not 
hi ghly contaminated by no n-fossil galaxies. We refer reader 
to IRicotti & Gnedir] J2005I) for more details on how the ob- 
servational sample is selected. 

2.3. Survival Probability 

In order to test the fossil scenario, we need first to measure 
the probability for a low-mass dark matter halo to remain rel- 
atively intact (i.e., not to evolve significantly) from z ~ 8 until 
the present time, which Ricotti & Gnedin (2005) called the 
"survival probability." We use the L20 simulation to measure 
the survival probabilities for dark matter halos. Since in the 
dissipationless simulation we cannot identify halos with lumi- 
nous galaxies, we can only measure the survival probability as 
a function of the maximum circular velocity. We will use the 
mapping between circular velocities and galaxy luminosities 
provided by the LI run (see Figure [3jl to calculate from V max 
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FIG. 3. — The luminosity - maximum circular velocity relation for the dwarf 
galaxies in the run LI at z = 8. The thin dotted line marks the resolution limit 
of our sample at V max = 13 km/s (see text for details). 
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FIG. 4. — The survival probability, PsWnax), of all dark matter halos that 
end up within 100 kpc (solid line), 300 kpc (dashed line), and 1 Mpc (dotted 
line) from the parent halo in the L20 run. The error bars are Poisson only. 
(Reminder: the maximum circular velocity V raax is measured at z = 8.) 



to luminosity. 

We measure the survival probability by tracking evolution 
of the individual halos fr om z = 8 to the presen t time 7 using 
the method developed by Kravtsov et al. (2004). To track the 
halos we use 156 simulation outputs equally spaced in expan- 
sion factor from z = 14 to z = 0. 

The tracks then can be used to estimate the fraction of halos 
existing at z = 8 which remain distinct, bound objects at z = 0. 
For detailed comparison with observations, it is useful also to 
differentiate halos with respect to the distance to the Milky 
Way-sized progenitor, as significant morphological segrega- 
tion is known to exist in the Local Group. Thus, we define 
the survival probability of halos within a given distance d (at 
z = 0) from the parent halo as 



iW,V ma x) = 



jy H (8^o,y max ) 

/VH(z=8,V max ) ' 



(1) 



where /Vh(8 — * 0, V max ) is the number of halos that exist both 
at z = and z = 8 within distance d (measured at z = 0) from 
the parent halo as a function of V max , and Nh(z = 8, Vmax) is the 
total number of halos at z = 8 located in the region of space 
(the "Lagrangian region") that ends up as a sphere of radius d 
centered on a parent halo at z = 0. 

It is important to stress that in equation Q the distance d 
is measured at z = 0, while the maximum circular velocity 
Vm ax is measured at z = 8. Operationally, in order to find all 
halos within the z = 8 Lagrangian region of sphere of radius 
d around the Milky Way-sized halo at z = 0, we tag all dark 
matter particles within the sphere and identify all dark matter 
halos at z = 8 located in the region of space occupied by the 
tagged dark matter particles. 

Note that in this study we calculate and apply the sur- 
vival probability to the small-mass dwarf objects of Vmax 
30 km/s. Such galaxies are dark matter dominated and effects 
of baryon dissipation on the inner mass distribution, V m ax, and 
on the survival of halos at small distances to the host (e.g., 
iMaccio et all2005l) should be small. 

As an example, Figure |4] shows the survival probability Ps 
as a function of V ma x(z = 8) for all dark matter halos from the 
L20 run. The survival probability ranges from a few to w 25 
per cent and there is a clear trend with the distance to the 
parent halo. Our measurement of Ps(V max ) becomes highly 
noisy for V ma x > 35 km/s due to small numbers of halos. 

7 In the next section we discuss the sensitivity of our results to the specific 
choice of the redshift at which fossils are identified. 



2.4. Defining a Fossil 

The survival probability is of course defined for all dark 
matter halos, not only for those that host fossil galaxies. In 
order to compute the survival probability for fossils only, we 
need an operational definition of a reionization fossil appli- 
cable to dissipationless simulations. We consider a dwarf 
galaxy to be a fossil if it experiences little or no star for- 
mation after reionization. This would be the case if the 
galaxy does not accrete fresh gas to fuel its continuing growth 
and star formation. The accretion of gas on dark mat - 
ter halos is controlled by the filtering mass (Gnedin 2000), 
which is easily calculable for a representative volume of 
the universe, but is not well defined for a region that col- 
lapses into a Milky Way-sized galaxy. As we noted above, 
such region can be reionized considerably earlier. Never- 
theless, after the region is ionized, the characteristic mass 
and circular velocity below which accretion is suppressed 
should be in the ra nge of ~ 30-50 km/s, a s was shown 
by many studies llThoul& Weinberg! II 996t lOuinnetalJ 
| 1996llNavarro & Steinmetzll997tlKitavama & rkeuchil2000t 
Gnedin 2000; Diikstra et al. 2004). Therefore, we define a 
reionization fossil as a halo that never increased its maximum 
circular velocity above a given fixed threshold Vf (note that 
we use Vf here to distinguish the threshold circular velocity 
for fossils from the filtering velocity, Vf), 



max[V ma x(z)] < Vrf. 



(2) 



This definition can be easily implemented for the dissipation- 
less simulations we use, as we track the evolution of V ma x for 
all halos. 

Figure [5] shows the dependence of the survival probability 
on this parameter for halos with V max = 20 ± 2 km/s. In the 
limit of large V r f, for example, 14% of all halos present at z = 8 
survive as separate object by z = for the d < 1 Mpc sample, 
while the rest grow into larger objects or get disrupted before 
z = 0. For objects closer to the parent halo the survival proba- 
bility decreases, since merging and disruption are expected to 
be more efficient in the higher density environment. 

The survival probability levels off for Vf > 50 km/s, be- 
cause if the value of the threshold V r f is set too high above 
the value of V ma x f° r the halos (20 km/s for Fig.|5Jl, it becomes 
irrelevant. Note also that the survival probability is a weak 
function of the filtering velocity for Vf > 30 km/s. Our results 
thus are not particularly sensitive to the specific choice of Vf 
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FIG. 5. — The survival probability Ps as a function of the threshold pa- 
rameter V r f (gray (me) for halos with V max = 20±2km/s from the L20 run. 
Three lines correspond to the three halo subsamples: halos within 1 Mpc (dot- 
ted line), 300 kpc (dashed line), and 100 kpc of the parent halo (solid line), 
respectively. 



as long as it is in this range. We demonstrate this explicitly 
below (see Fig.0. In what follows, we adopt Vrf = 30km/s as 
our fiducial value, as it is sufficiently small to ensure that ha- 
los that we identify as fossils have minimal accretion of fresh 
gas after reionization. 

Figure |6] shows the survival probability as a function of 
V max (z = 8) (similar to Fig. [4} but for the reionization fossils 
defined using equation with Vrf = 30 km/s. Note that Ps ap- 
proaches zero for V max approaching the value of the threshold 
Vrf = 30 km/s, simply because the definition does not allow ex- 
istence of such fossils. At V max < 15 km/s the incompleteness 
of the d < 100 kpc sample becomes non-negligible, and the 
survival probability at lower values of the maximum circular 
velocity becomes uncertain by up to a factor of two. This un- 
certainty does not significantly affect our results since most of 
our galaxy sample has V max > 15 km/s. 

The survival probability ranges from a few per cent for ob- 
jects located close to the Milky Way host to ~ 15% for more 
distant halos. Thus only a small fraction of objects at z = 8 
would survive to be true fossils at the present epoch. The sur- 
vival probability distribution shown in Figure|6]is the main in- 
gredient of our model for the luminosity function and spatial 
distribution of fossil galaxies presented in the next section. 

3. RESULTS 

Given the survival probability of fossils as a function of 
Vmax and d from the L20 run, we can use the galaxy forma- 
tion run LI to compute the luminosity function of fossils as a 
function of distance from the parent galaxy. Specifically, we 
use the distribution of high-redshift dwarfs in the Ly - V max 
plane shown in Fig.|3]to compute the cumulative luminosity 
function of fossils at z = within a given Lagrangian region 
with mass Mlr: 



JVlrO Lv) ■■ 



Mi 



LR 



1.2 x 1O U M 



1O 7 L 



dL 



^V max P s (V max ) 



d 2 N 8 (L,V m . dx ) 
dLdV max 



(3) 



where 1 .2 x 10 1 1 M Q is the mass of the computational volume 
of the LI run and d 2 Ns,(L,V max ) / dL/ dV max is the luminosity 
function of high redshift dwarfs at z = 8 in the L 1 run per unit 
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FIG. 6. — The same as in Fig.|4] but for the reionization fossils according to 
the definition J5J. The solid gray line shows the survival probability for the 
d < 100 kpc when the correction for incompleteness is taken into account. 



V max . For the small number of objects we have in the LI run, 

d 2 N s (L,V max ) 



dLdV m 



i 



5(V ma x-V max )<5(L v -4), 



where V max and L J V are the maximum circular velocity and the 
V-band luminosity of a dwarf galaxy /'. 

The upper limit of 10 L^, in the integral over the luminos- 
ity is motivated by Ricotti & Gnedinl (j2005), who found that 
no fossils exist at higher luminosities in the Local Group. The 
limits of the integral over V max are formally over all existing 
halos, but the survival probability is non-zero only over a fi- 
nite range of 1 3 km/s < V max < Vrf due to our definition of the 
reionization fossil and our completeness limit. 

In order to define the amplitude of the cumulative lumi- 
nosity function for the Local Group, we need to know the 
mass of the Lagrangian region that ends up within a given dis- 
tance from the parent galaxy (the Milky Way or Andromeda). 
To estimate this mass, we use the dynamical m odels of the 
Milky Way and Andromeda from Klvpin et al. (2002). For 
their preferred models (Ai or Bi for the Milky Way and Ci 
for Andromeda), we find the total mass with 100 kpc from 
the Milky Way or Andromeda (summed over these two galax- 
ies) to be 1.4 x 1O 12 M , and the total mass within 300 kpc 
to be 2.7 x 1O 12 M . For larger distances, extrapolation is 
required. In the L20 run, about 50% more mass is between 
300 kpc and 1 Mpc around the parent halo compared to the to- 
tal mass within 300 kpc. This number is also consistent with 
the simple extrapolation of NFW profiles around the Milky 
Way and Andromeda. Thus, if the two galaxies were iso- 
lated, the total mass of the Local Group within 1 Mpc would 
be » 4.0 x 10 12 M Q . However, since the distance between 
the Milky Way and Andromeda is about 780 kpc, there is a 
substantial overlap between the two spheres of 1 Mpc in ra- 
dius centered on the two galaxies, and most of the mass of 
both galaxies is within the overlap region. Although an ac- 
curate estimate is difficult, we adopt the mass 3.0 x 1O 12 M 
as a reasonable estimate of the total mass of the two galaxies 
within 1 Mpc from each of them. 

We first check the sensitivity of our results to adopted val- 
ues of free parameters: the starting redshift z s at which we 
identify the future fossils and the fossil threshold Vrf. Indeed, 
although z = 8 is a reasonable choice for characteristic pre- 
reionization epoch, there is nothing special about it. This spe- 
cific value is simply the epoch when the LI run was stopped in 
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FIG. 7. — The cumulative luminosity functions of fossil dwarfs located 
within 300 kpc from the parent halo at z = as a function of the starting 
redshift (top panel) and the fossil threshold parameter Va (bottom panel). 



iRicotti & Gnedinl J2005I) for physical and computational rea- 
sons. In addition, depending on the particular evolutionary 
histories of the real Milky Way and Andromeda, the Local 
Group region could have been reionized at different redshifts. 

The top panel of Figure^ shows the cumulative luminosity 
functions for fossils located within 300 kpc of the parent halo 
at z = 0, calculated using survival probabilities derived from 
the L20 run. Four lines correspond to the four values of start- 
ing redshift: z s = 7, 8, 9, and 10. Namely, the solid line in 
Fig. is obtained using the equations above, while the dot- 
ted line (z = 7) is obtained from the same equations but with 
z = 8 replaced by z = 7 everywhere. The figure shows that the 
dependence on the specific value for the starting redshift for 
measuring the survival probabilities is quite weak for z s < 9. 
In making this figure, we used the same Ly - V mdx distribution 
from Fig.|3j since we cannot rerun the LI simulation for dif- 
ferent reionization redshifts (due to its expense). However, a 
preliminary analysis of different galaxy formation simulation 
(Tassis et al. 2006, in preparation) indicates that redshift de- 
pendence of the L - V max relation for high-z dwarfs is given 
byia Vmax(l +z) -1 ' 5 - Depending on how different the red- 
shift of reionization is from our fiducial value of z = 8, the 
difference in the relation can be approximated by this redshift 
scaling. Thus, for example, if we change redshift from z = 8 
to z = 10 the difference in luminosity for a given V max is only 
35%. If however, the change is to z = 15 the difference is a 
factor of two. This uncertainty can therefore be assumed as 
uncertainty of a factor of two in the luminosity we assign to 
the fossil dwarfs. 

At higher redshifts the dependence is much stronger. We 
can of course expect that as we increase the starting redshift 
to z s ~ 10-20 the number of halos that do not change appre- 
ciably to the present epoch becomes very small. Our choice 
for the lower value of z s is fairly straightforward. The uni- 
verse is known to be fully reionized at z ~ 6, so z = 7 is quite 
close to the lowest possible value. Indeed, it is highly unlikely 
that the vicinity of the Local Group should be one of the last 
patches of neutral gas in the universe to be re-ionized at z < 7. 
The upper limit of z ~ 9 is less constrained, but is generally 
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FIG. 8. — The cumulative luminosity functions of fossil dwarfs located 
within 1 Mpc (top panel), 300 kpc (middle panel), and 100 kpc (bottom panel) 
from the parent galaxy at z = 0. The black lines with error-bars are the ob- 
servational data for the Local Group dwarfs, while the light gray bands show 
model predictions with la error-bars, uncorrected for incompleteness. The 
dark gray line in the bottom panel shows the model prediction corrected for 
the incompleteness. 



a reasonable number to expect a region of about l/i" 1 Mpc 
across to be reionized by the UV radiation from a Milky Way 
progenitor. Of course, other sources of radiation may cause 
even earlier partial or complete ionization of the Lagrangian 
region around the Milky Way or Andromeda. The predicted 
abundance of fossils would be correspondingly smaller, as we 
discuss in the next section. The bottom panel of Figure [7] 
shows dependence of the luminosity function of fossils on the 
adopted value of Va . The dependence is weak and the ampli- 
tude of the luminosity function changes by less than a factor 
of two for Vrf from 25 km/s to 35 km/s. 

Figure [3] presents the main result of this paper: the cumu- 
lative luminosity function of reionization fossils in the Local 
Group within three different distances from the parent galaxy. 
The o bservational data for the Local Group dwarfs ar e taken 
from lMate'ol Jl9 98) and McConnachie & IrwTnl J2005I) . as de- 
scribed in § 12.21 Since observational uncertainties on the dis- 
tances (from the Milky Way) to some of Andromeda satel- 
lites are comparable to their distance to Andromeda, there 
is an uncertainty on whether a given dwarf galaxy is within, 
say, 100 kpc from Andromeda. To take this uncertainty into 
account, we construct 100 random realizations of the Local 
Group in which dwarf spheroidal galaxies are distributed over 
the distance (from the Milky Way) with a log-normal distribu- 
tion of the width corresponding to the distance measurement 
error. 

The errorbars that we assign to the observed luminosity 
function include both the Poisson fluctuations and the error 
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FIG. 9. — The radial distribution of reionization fossils with Ly > 10 6 Lq 
(top panel) and > 10 s Lq (bottom panel). The light gray bands and dark 
gray lines show theoretical expectations without and with correction for in- 
completeness. Black lines with error-bars show the observation al data on 
Local Group dwarfs identified as reionization fossils by Ricotti & Gnedin 
1 2005), while the black dashed line show all (whether a fossil or not) Lo- 
cal Group dwarfs. 



measured from the dispersion between 100 random realiza- 
tions of galaxy distances. The predictions of our model for 
the fossil galaxies is shown by the hatched band and includes 
Poisson errors in the number of galaxies from the LI run and 
Poisson errors in the survival probability due to the limited 
number of halos in the L20 run. Since we have only one dis- 
sipationless simulation, additional uncertainty due to cosmic 
variance should exist as well. However, we expect it to be 
smaller or comparable to the Poisson errorbars shown here. 
The dark gray line in the bottom panel of Fig. [8] shows the 
mean value of the theoretical prediction with the incomplete- 
ness correction included. 

We can draw two straightforward conclusions from Fig- 
ure|8] 

1. The expected number of fossils with Ly > 10 6 Lq is 
in reasonable agreement with the observed luminosity 
function of the Local Group dwarfs, which are likely to 
be fossils, for all distances from the parent galaxy. 

2. The predicted number of fossils with 10 5 L Q < Ly < 
10 6 L Q is also in reasonable agreement with observa- 
tions within 100 kpc, but is larger by a factor 2 to 4 for 
d < 300 kpc and d < 1 Mpc. 

These conclusions are not sensitive to our incompleteness cor- 
rection, as it is smaller than our estimated error bars. 

Figure |9] compares the predicted radial distribution of the 
fossil galaxies for two luminosity ranges (Ly > 10 6 Lq and 
Ly > 10 5 L Q ) to the distribution of observed Local Group 
dwarfs. It is interesting to n ote that the distr i bution of ob- 
served galaxies identified by Ricotti & Gnedin (2005) as fos- 
sils is quite similar to the radial distribution of all dwarf galax- 
ies. There is thus no particular spatial bias exhibited by such 
galaxies. It is not clear whether a strong spatial bias should 
exist. On the one hand, the early forming high-redshift galax- 
ies should be spatially biased with respect to their host. On 



the other hand, the galaxies located close to the host would be 
accreted early and thus may have a lower chance to survive 
until the present. 

The figure shows that the predicted radial distribution for 
the brighter fossils (Ly > 10 6 Lq) agrees with observations. 
For lower luminosities, however, the predicted and observed 
radial distributions are quite different. Although the number 
of predicted fossils at d < 100 kpc is similar to the observed 
abundance, the observed dwarfs show a somewhat more ra- 
dially concentrated distribution. The difference is not large 
however and is not statistically significant. At larger dis- 
tances, on the other hand, the model predicts far more (a fac- 
tor of 3-4) fossils with Ly > 1O 5 L0 than is observed. There is 
thus a marked deficit of the observed galaxies at the faintest 
luminosities at distances larger than 100 kpc. These conclu- 
sions are the same if we consider radial distribution for An- 
dromeda or the Milky Way only, although the statistics is 
poorer in this case. We discuss possible explanations for the 
discrepancy in the next section. 

4. DISCUSSION AND CONCLUSIONS 

In this paper we estimate the expected abundance and spa- 
tial distribution of the reionization fossils: dwarf galaxies that 
formed most of their stars around z ~ 8 and evolved only little 
since then. We show that ~ 5 - 15% of the objects existing 
at z ~ 8 do indeed survive until the present in the Milky Way 
like environment without significant evolution. Because such 
galaxies form their stellar systems early during the period of 
active merging and accretion, they should have spheroidal 
morphology and their observed counterparts are most likely 
to be identified among the dwarf spheroidal galaxies. The 
spheroidal morphology of the fossil galaxies is thus not re- 
lated to the tidal heating by the nearby massive galaxies and 
they can be found even at distances as large as a megaparsec 
from their host. This could possibly explain existence of dis- 
tant dwarf spheroidal galaxies in the Local Group, such as 
Cetus, Tucana, and KKR25. 

We show that both the expected luminosity function and 
spatial distribution of dark matter halos which are likely to 
host fossil galaxies agrees reasonably well with the observed 
distribution of the luminous (Ly > 10 6 L Q ) Local Group fossil 
candidates near the host (d < 200 kpc). The predicted abun- 
dance is substantially larger (by a factor of 2-3) for fainter 
galaxies (Ly < 1O 6 L ) at larger distances (d > 300 kpc). This 
discrepancy can have at least three plausible explanations. We 
list them and the possible tests below. 

Observations are wrong. Observational incompleteness for 
faint, Ly < 10 6 L©, fossils may potentially account for 
the difference we fin d, if it can be as high as a factor 
of 3 at d > 300 kpc. Willma n et all JJ004) found that 
a factor of 3 incompleteness is possible, although that 
number was on the very edge of their estimates. The 
same incompleteness factor would have to apply to the 
distribution of Andromeda satellites, which seems to be 
less likely. Future observational searches (for example, 
from SDSS) for dwarf spheroidals between 200 kpc and 
1 Mpc from the Milky Way or Andromeda, or observa- 
tions of nearby galaxy groups with the next generation 
large telescopes should be able to test this possibility. 

Theory is wrong. The abundance of fossils galaxies may be 
overestimated by our model. This is possible for ex- 
ample if reionization occurs much earlier than we as- 
sumed: for example, at z > 10, as may be suggested 
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by the first year WMAP data llKogut et al]l2003l) . In 
this case, the number of small-mass objects that are 
able to form stars before reionization should be much 
small er. However, if most of dwarf galaxies that 
Ricotti & Gnedin (2005) identified as reionization fos- 
sils are confirmed to be such by future observations, 
then this work will places a constraint of z < 10 on the 
reionization history of the region of the universe that 
collapses into the Local Group. 

Spatial bias. If we take the observational data as complete 
and at face value, they indicate that lowest luminosity 
fossils (Ly < 10 6 L Q ) are located preferentially near the 
parent galaxy. At first glance, such a conclusion seems 
counterintuitive, since the feedback (both radiative and 
kinetic) from the progenitor of the main galaxy is ex- 
pected to suppress star formation in nearby galaxies. 
However, galaxies are also biased, more so at high red- 
shift, so dwarf galaxies located closer to the progen- 
itor of a large spiral galaxy will start forming sooner 
than similar galaxies in more distant regions of space. 
Which of the two opposing effects (feedback vs bias) 
wins is unclear, but should be tackled with the future 
ultra high-resolution cosmological simulations. 

The results of this study indicate that existence of the fossil 
galaxies in the Local Group, in which most of the stars formed 
during the earliest stages of galaxy formation, is plausible. 



The abundance of fossils predicted by our model would con- 
stitute as much as 20-40% of all the luminous Local Group 
dwarfs. The improvement in resolution of cosmological sim- 
ulations should improve our estimates. Observationally, iden- 
tification and studies of the fossil galaxies can provide an ex- 
citing opportunity of studying the properties of star formation 
in primeval galaxies by direct observations in our own cosmo- 
logical backyard. 
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